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ABSTRACT 

Context. Multi-wavelength galaxy number counts provide clues on the nature of galaxy evolution. The interpretation 
per galaxy type of the mid-IR faint counts obtained with ISO and Spitzer, consistent with the analysis of deep UV- 
optical-nearlR source counts, provide new constraints on the dust and stellar emission. Discovering the nature of new 
populations, as ultra- luminous (>1O^^L0) infrared galaxies (ULIRGs), is also crucial for understanding galaxy evolution 
at high redshifts. 

Aims. We first present the faint galaxy counts at 12 /im from the catalogue of the KO-ESO-Sculptor Survey {ISO -ESS) 
published in a companion article (Seymour et al. 2007), which go down to 0.24 mjy after corrections for incompleteness. 
We verify the consistency with the existing ISO number counts at 15 fim. Then we analyse them as well as the 24 fim 
{Spitzer) faint counts to constrain the nature of ULIRGs, the cosmic star formation history and the time-scales for 
mass build-up. 

Methods. We show that the "normal" scenarios in our evolutionary code PEGASE, which had previously fitted the deep 
UV-optical-nearlR counts, are unsuccessful at 12 and 24 /xm. We then propose a new scenario of ULIRG adjusted 
to the observed cumulative and differential 12 /im and 24 counts and based on observed 12 and 25 fim IRAS 
luminosity functions and mid-IR colors from PEGASE. 

Results. We succeed in modelling the typical excess simultaneously observed at 12 /im, 15 {ISO) and 24 fim {Spitzer) 
in the cumulative and differential counts by only changing 9% of normal galaxies (1/3 of the ellipticals) into ultra- 
bright dusty galaxies evolving according to the scenario of ellipticals, and interpreted as ULIRGs. These objects present 
similarities with the population of radio-galaxy hosts at high z. No number density evolution is included in our models 
even if rare occasional starbursts due to galaxy interactions remain compatible with our results. 

Conclusions. Higher spectral and spatial resolution in the mid-IR together with submmillimeter observations using the 
future Herschel observatory will be useful to confirm these results. 

Key words. Cosmology: surveys; Galaxies: evolution - spiral - infrared - photometry 



1. Introduction 

The mid-infrared extra-galactic source counts provide clues on the evolution of galaxies at high z and allow us to follow 
the cosmic star formation history up to when the Universe was only one quarter of its present age. The high sensitivity of 
the infrared deep surveys observed with ESA's Infrared Space Observatory ISO (Kessler et al. 2003), and more recently 
with Spitzer Space Telescope (Werner et al. 2004) offers an unique opportunity to study the obscured star formation 
process through the emission of grains heated by young stars and possibly by an active nucleus, thus reaching galaxies at 
their earliest epochs. Interactions of galaxies are known to explain the huge IR emission detected in Ultra Luminous Infra 
Red Galaxies (ULIRG; Soifer et al. 1984). One crucial issue is whether the ULIRGs provide observational evidence that 
galaxy interactions play a fundamental role in galaxy evolution. Moreover, because the infrared luminosity depends on 
the dust mass accumulated from stellar ejecta, it provides complementary diagnostics of past star formation which may in 
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turn be used to model the star formation history. Over the years, a long series of deep surveys from the UV to the optical, 
down to the extreme depth of the Hubble Deep Field North (HDF-N) at B = 29 (Williams et al. 1996), have constrained 
the direct stellar emission in terms of cosmology and scenarios of star formation. Fioc and Rocca-Volmerange (1999a) 
derived a set of galaxy populations fitting the multi- wavelength (UV-optical-ncar infrared) deep galaxy counts dominated 
by stellar emission. This set defines the evolution scenarios of eight various galaxy types and their number fractions. In 
the mid-infrared, galaxy light is dominated by the dust emission from grains in the form of graphite, silicates or Polycyclic 
Aromatic Hydrocarbons (PAH; Puget & Leger 1989). Because different time-scales characterise the respective star and 
dust emissions, the results of optical and mid-infrared source counts may significantly differ. 

However two difficiiltics hamper the interpretation of data. One is the lack of homogeneity between the various surveys: 
due to large-scale galaxy clustering, the statistical properties derived from deep pencil beam surveys suffer from "cosmic 
variance" (Somerville et al. 2004), and might thus differ from the analyses on large area surveys. Another difficulty is 
the variety of sources (star bursts due to mergers, normal evolved galaxies, AGNs) and their intrinsic evolution. As an 
example, the time-scale of star formation associated to galaxy interactions is significantly shorter (^lO^yrs) than the 
star formation time-scales of galaxy populations (> lO^yrs, depending on spectral type) observed in deep surveys. We 
consider the contribution of AGNs to be minimal at the flux density range explored here, but the impact of an embedded 
hidden AGN is discussed below. 

Our large area ISOCAM ESO-Sculptor Survey (ISO-ESS) at 12 /xm, published in the companion paper (Seymour et 
al, 2007), here tackles these various limitations by covering a significant area, and by using the new code PEGASE.3 
(Fioc et al. 2007) which coherently predicts the evolving stellar and grain emissions from evolved galaxies as well as young 
starbursts. It is able to predict both starlight and dust emission from the optical to the thermal infrared, by taking into 
account the transfer and the reprocessing of light in the different wavelength domains. A large variety of star formation 
time-scales is considered in our evolutionary templates but no evolution of the number density of galaxies is included in 
the model. 

Several other major surveys in the mid-infrared have been performed with the ISOCAM camera which also provide 
deep galaxy counts. The largest survey is ELAIS (Rowan- Robinson et al. 1999, 2004), which covers 12 sq. dcg. in the fiux 
range 0.45-150 mJy at wavelengths of 6.7 /zm (LW2) and 15 /zm (LW3); the corresponding galaxy counts were published 
by Serjeant et al. (2000). Following the preliminary surveys of Taniguchi et al. (1997) and Oliver et al. (1997), several 
major surveys have also been performed which provide deep galaxy counts. Among them arc the ISO 15 /im observations 
of the Lockman Deep Field and the Marano-ROSAT Ultra Deep Field (Aussel et al. 1999, Elbaz et al. 1999) and the 
Lockman Shallow Field (Flores et al. 1999). Even deeper surveys have been obtained in the LW2 and LW3 filters, centered 
on Abell cluster 2390 (Altieri et al. 1999), the HDF field, and others fields covering areas of 2.5 arcmin in radius (Oliver 
et al. 2000, 2002), and of 16 arcmin^ (Sato et al. 2003). These various surveys yield galaxy number counts which are 
in reasonable agreement, and detect an excess in the number of galaxies at faint fiuxes (below 1 mJy) which is often 
interpreted as an increase in the star formation rate with look-back time. For example, Pozzi et al. (2004) fit the 15 /xm 
number counts by introducing a population of evolving starbursts, based on the local star-forming prototype M82 (Silva 
et al. 1998), which undergoes very strong luminosity or density evolution parameterized as (1 -|- z)", with a ~3-4. 

All the differential 24 /xm counts also show a systematic departure from an Euclidean Universe at ffuxes fainter than 
^ 0.5 mJy (Rodighiero et al. 2006), with a peak at ^ 0.3 mJy, similarly to the 15 /im and 12/im observations. Below 
60/xJy, the 24 /xm galaxy number counts obtained by the Spitzer/MlPS deep surveys (Marleau et al. 2004, Papovich 
et al. 2004, Chary et al. 2005) face the problem of confusion limit due to extragalactic sources. Several interpretations 
propose huge evolution factors in luminosity and/or density. In the Chandra Deep Field South survey, the sample of 2600 
Spitzer/MlPS sources brighter than 80/iJy is analyzed with the comoving IR energy varying as (1 -|- z)^-^^^-'^ (Le Floc'h 
ct al. 2005). The origin of such a high evolution is however not described. More recently, Caputi et al. (2006) analyzed 
the stellar populations of the Spitzer/MlPS 24 /um galaxies in the GOODS/CDFS from Kg images. They show evidence 
for a bump in the redshift distribution at 2; = 1.9, induced by a significant population of galaxies with PAH emission. 

The combined ISO-ESS presented here has comparable depth (^ 80% completeness at ^ 0.7 mJy) and surface area 
(680 arcmin^) at 12 /im as the intermediate surveys (Lockman Deep and Marano Deep fields) performed at 15 /im as part 
of the ISOCAM Guaranteed Time Extragalactic Surveys. The specificity of the ISO-ESS survey is to provide deep galaxy 
counts at ~ 12 ± 3.5 fim (LWIO filter), in a wavelength range where the spectral energy distribution (SED) is dominated 
by the signatures of PAH emission; it thus strongly constrains the evolution process. The target field is located within 
the ESO-Sculptor Survey (de Lapparent et al. 2003, 2004), for which deep optical BVRc magnitudes up to Rc < 23.5 
and spectroscopic redshifts at Rc < 21.5 have been obtained (Arnouts et al. 1997, Bellanger et al. 1995a). 

We then model the 12 /xm number counts using the new version PEGASE.3 of the "Projet d'Etude des GAlaxies 

par Synthcsc Evolutive" (Fioc et al. 2007; Fioc & Rocca-Volmerange 1997, 1999b; sec also www2.iap.fr/pegase), which 
coherently complements the UV-optical-NIR emission from stars and gas with the mid- and far-infrared emission from 
dust. The specific goal of our analysis is to predict the mid- infrared number counts using firstly the same scenarios of 
galaxy evolution by type, with the same number fractions and the same total number densities of galaxies as used for 
the successful predictions of the UV-optical-NIR deep counts by Fioc & Rocca-Volmerange (1999a). Secondly, if needed, 
other scenarios arc proposed to model the population of ULIRGs able to reproduce the systematic departure observed 
around 0.3mJy in mid-infrared surveys. 

Section 2 describes the parameters of the ISO-ESS field observed with the large pass-band ISOCAM LWlO/12/im 
filter and summarises the data analysis required to extract a catalogue of sources presented in a companion article 
(Seymour et al. 2007). The corresponding cumulative and differential galaxy counts at 12 /xm are presented in section 3. 
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We resume in section 4 the parameters of two models of evolution scenarios and the observed IRAS luminosity function 
at 12 /xm, tentatively used to model galaxy counts. The respective fits of the cumulative and differential counts of the 
ESO-ESS survey by the two models are compared in section 5. The best model 2 identifies a population of ultra bright 
ellipticals interpreted as distant ULIRGs. Applied to 24 /im number coimts from the deep Spitzer/MTPS surveys, we show 
in Section 6 that the same population of ultra bright ellipticals also reproduces the 24 /xm differential excess, making 
robust our model 2. Section 7 presents the cosmic star formation history resulting from the best fit, taking into account 
the respective contributions per galaxy type. Section 8 discusses the stellar and dust masses of the revealed ULIRG 
population, and the possibility of an hidden AGN. The final section presents our conclusions. 

2. Observations and data reduction 

2.1. The ESO-Sculptor Survey field 

The selected field, the ESO-Sculptor Survey (ESS) of faint galaxies, is described in its complete version by de Lapparent 
et al. (2003). Located near the Southern Galactic Pole, the observations for the ESS were performed as an ESQ key- 
programme, thanks to guaranteed time on the ESO NTT and 3.6m telescope. Deep CCD Johnson B, V and Cousins 
i?c magnitudes for nearly 13000 galaxies to ~ 24 were obtained over a continuous area of ~ 0.37 deg (Arnouts et al. 
1997). Multi-object spectroscopy has also provided redshifts and flux calibrated spectra over a sub-area of ~ 0.25deg^ for 
617 galaxies with < 20.5 (92% complete) and 870 galaxies with Rc < 21.5 (52% complete). The optical star/galaxy 
separation was performed using the "stcllarity" index from the SExtractor software (Bcrtin & Arnouts 1996). The 
optical spectra were classified using a Priiicnpal Component Analysis (Galaz & do Lapparent 1998). The optical luminosity 
functions were then measured per galaxy spectral type (de Lapparent ct al. 2003), and lead to the detection of a marked 
evolution in the spiral galaxy populations, characterised by an excess of late spiral and irregular galaxies at z ~ 0.5 — 1.0 
(de Lapparent et al. 2004). 

2.2. The ISO survey on the ESS field (ISO-ESS) 

The ISO observations of the ESS field have been performed with the raster mode CAMOl of the ISOCAM camera 
(Ccsarsky ct al. 1996) on board of the Infrared Space Observatory ISO (Kessler ct al. 2003). The broad-band LWIO filter, 
with reference wavelength Ard = 12 /im and covering the 8.5 15.5 /im interval, has been built for a direct comparison 
with the 12 /xm IRAS filter (Moneti et al. 1997). The pixel field of view (PFOV) with ISOCAM is a 6 arcsec square and 
the adopted integration time was 5.04 sec. per exposure. Ten raster maps were built within the total on-target time of 
14h. All the rasters have M x = 8 x 8 pointings each offset by dM = dN = 60 arcsec along the axis of the detector. 
The number of exposures per pointing is A'exp = 13 and the number of stabilization exposures is iVgtab = 10. The total 
area of the survey is ~ 680 arcmin^ with a maximum exposure time of 1200 sec. and an average exposure time of ~ 660 
sec. per sky pointing. 

The ISO target field has been selected in the region of the ESO-Sculptor survey where the cirrus emission is minimal 
(see companion article). The overlap between the 2 surveys is large, as it represents 90% of the ISOCAM survey area, 
and 75% of the ESO-Sculptor spectroscopic area. 

2.3. Summary of the data processing 

This section is a summary of the data analysis presented in the companion article (Seymour et al. 2007). The various 
steps of data processing are: (1) deglitching (i.e. cosmic ray subtraction), (2) correction of the long transient behavior 
etc. performed by the PRETI software of Stark et al. (1999), and (3) source detection above a noise map, performed 
using a wavelet analysis technique. The detailed adaptation of the PRETI method to the ISOCAM data analysis was 
published by Aussel et al. (1999). We fine-tuned the astrometry of the 12 /xm sources by comparison with matched 
objects in the 2MASS catalogue, which led to a maximum change of 0.3 arcsec in the position of the 12 /xm sources. 
The subsequent r.m.s. offset between the ISOCAM 12 /im and the ESS coordinates is 2 arcsec, with no systematic offset. 
The star/galaxy separation was performed by a detailed analysis of the near-infrared (NIR) and optical colour-colour 
diagrams. We then derived an empirical flux calibration from a subset of stars in our sample. This calibration was based 
on fitting the optical-NIR counterparts of the stellar data with stellar spectrum templates compiled in the PEGASE 
library and then using the optical/infrared relations from IRAS data to predict the true infrared fluxes of stars in our 
field. This procedure took advantage of the deep photometric ESS survey in the optical. Our flux calibration is robust as 
it was carried out using two different colour-colour relations and empirical calibration using well-known stellar templates. 
The final catalogue of 142 sources with ISO fluxes is listed in the companion article , and contains 22 stars and 120 
galaxies. As expected, only a small fraction of the MIR sources are stars. 

2.4. Large-scale distribution of the ISO-ESS sources 

With an average exposure time of 11 minutes by pointing, our ISO-ESS survey is as deep as the Lockman Hole Deep 

survey from the /S'OCVIM Guaranteed Time Extragalactic Survey (Elbaz ct al. 1999), but over a 33% larger area of ~ 680 
arcmin-^. Despite its significant area, the ISO-ESS map reveals the inhomogeneity of the projected distribution of MIR 
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Fig. 1. The distribution by types of the (R<20.5) ISO-ESS galaxies, derived from the optical spectral classification (de 
Lapparent et al. 2003). The dominant population is Sb-Sbc spirals, confirming that ellipticals are rarely detected at 
12 ^m. Note however that among the 104 galaxies detected at 12 /j,m in the area common to both surveys, only 77 have 
ESS spectroscopic information (hence redshift and spectral type) , leaving place for the discovery of new populations and 
signatures of evolution. 



sources. The large-scale distribution of galaxies in the ESS shows the alternation of sharp walls and voids (Bellanger & 
de Lapparent 1995b). For Ho=65 km s~^ Mpc~^, de Lapparent & Slezak (2007) measure a comoving correlation length 
of ~ 5.5 Mpc at the median redshift z ^ 0.3 of the ISO-ESS galaxies; the right ascension transverse extent of the survey 
(~ 11 Mpc) thus corresponds to approximately twice the galaxy correlation length. Because the MIR sources are likely 
to follow the clustering of optical galaxies, large-scale fluctuations are naturally expected. These are however smoothed 
out when calculating the ISO-ESS galaxy counts summed over the whole redshift range of the survey. 



3. Faint galaxy counts at 12 fim 

The ISO-ESS sample consists of sources with a detection level equivalent to 5a and is complete to 1 mJy. At fainter 
flux densities, the correction for incompleteness in the interval 0.24-1.0 mJy is computed by two independent methods 
respectively based on stars and galaxies in the optical. Both approaches take advantage of the deep ESS optical survey 
which contains counterparts to all sources detected by ISOCAM down to 0.24 mJy (in the common area to both surveys), 
and yield incompleteness corrections which are in good agreement (see Fig. 7 in companion paper). Here, we correct the 
galaxy number counts using the incompleteness correction derived from the stars, as it is least affected by the large-scale 
clustering in the ESS sample (see sect. 2.4). Fig.[T]shows the distribution by types of the 77 /5'0-ESS galaxies which have 
measured redshifts from the optical survey (de Lapparent et al. 2003). It shows that spirals Sb and Sbc are the most 
numerous galaxies detected in the MIR. It also confirms that gas-poor normal ellipticals are essentially undetected. The 
ESS spectral classification used the code PEGASE.2 (Fioc & Rocca-Volmerange 1997). However, the sample on Fig. [T] 
is incomplete, with only 77 galaxies having a measured redshift. Among the 104 ESS galaxies detected at 12 /im, the 
27 objects with no redshift measurement could be at higher redshift or belong to a new type. We therefore adopt the 
"optical" density fractions derived from the faint count analysis in the UV-optical-nearlR (~ 27%, 30%, and 43% for 
early-, intermediate-, and late-type galaxies respectively) by Fioc et al. (1999a), confirmed by de Lapparent et al. (2003) 
in the Sculptor field. 



3.1. Cumulative galaxy number counts N(> S^) at 12 ^m 

We then derive the number of detected ISO-ESS galaxies as a function of 12 /im flux density. The detections are binned 
in 11 flux density intervals, chosen so that each bin contains 10 — 12 more sources than the previous one. The resulting 
cumulative source counts N{> S^) (commonly referred to "integrated counts" in the literature), where Si, is the 12 /tm 
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Fig. 2. Cumulative 12 fim galaxy counts observed with the LWIO filter from the ISO-ESS survey, plotted as full squares 
with error bars. The other ISOCAM 15 /im surveys are also plotted: the compilation by Elbaz et al. (1999) as a solid line 
contour, which covers the Lockman Hole (Rodighiero et al. 2004), the Marano field (Aussel et al. 1999), and the HDF 
North and South fields (Altieri et al. 1999) to various depths; the number-counts from the ELAIS survey (La Franca et 
al. 2004) are shown as a dotted line. The 15 /im counts are corrected to 12 /im, see text for details. 



flux in mJy, derived after correction for incompleteness, are shown in Fig. [2]as full squares. The vertical error bars in the 
corrected cumulative counts are estimated by assuming Poisson fluctuations in the number of detected galaxies, and by 
taking into account the uncertainty in the incompleteness correction; the horizontal error bars are not plotted as they are 
smaller than the size of the symbols. The large error bars in the high flux point at Si, ~ 1.4 mJy result from the small 
number of galaxies (10) detected at these high flux densities: this is evidently due to the limited area of the ISO-ESS 
field. 

We overlay in Fig. [2] the results from the compilation of ISOCAM 15 /im surveys obtained in the LW3 filter as 
published by Elbaz et al. (1999; solid line): they cover the Lockman Hole, Marano and HDF North and South fields 
(Rodighiero et al. 2004; Aussel et al. 1999; Altieri et al. 1999) at various depths. We also plot for comparison the results 
from the ELAIS survey, also observed at 15 /im (La Franca et al. 2004, dotted line). When plotting the 15 /im number 
counts, the flux density is converted into the LWIO 12 /im band using the respective central wavelengths of the two fllters 
and assuming a flat spectrum in flux (that is, a constant product of the frequency by the flux density); the correction 
is small, a factor 12/15 = 0.8. We do not correct the counts from 15 /im to 12 /xm because the filters are wide and very 
close to each other in wavelength (they actually slightly overlap). 

Figure [2] shows that our data are in good agreement with the surveys compiled by Elbaz et al. (1999), including 
the deepest samples. In particular, the faint-end slope of the ISO-ESS cumulative counts is similar to that of the ultra 
deep survey obtained on the cluster-lens A2390 by Altieri et al. (1999) in their common flux density interval, despite 
the narrow pencil-beam geometry of the latter survey. The ISO-ESS counts are also consistent with the brighter galaxy 
survey from La Franca et al. (2004). 
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Fig. 3. The ISO-ESS 12 /im differential counts normalized to the Euclidean case, plotted as filled squares (note that 
-S'^-^dN/dS is positive). The other following 15 surveys are also shown: ELAIS-Sl (Pozzi et al. 2003, triangles); 
through the cluster- lens A239G (Altieri et al. 1999, 6-pointed open star), HDF North (Aussel et al. 1999, crosses); HDP 
South (asterisks), Marano PIRBACK Ultra Deep field (open squares), Marano Ultra Deep ROSAT field (points), Marano 
Deep field (open circles), all published in Elbaz et al. (1999); Lockman Deep field (4-pointed open star), Lockman Shallow 
field (3-pointed heavy cross), both from Rodighiero et al. (2004). See text for the details of the correction from 15 /xm to 
12 /xm. 



3.2. The Euclidean-normalized differential counts 

The differential galaxy counts -S^'^dN/dS{>0 ) of the ISO-ESS survey at 12 /im, normalized to those for an Euclidean 
Universe are presented as filled squares in Pig. [31 A wider binning than in Pig. [2] is adopted, with 23 — 25 galaxies per bin. 
The vertical error bars take into account the Poisson error in the number of actually detected sources, combined with the 
uncertainties in the incompleteness correction and in the fiux density of each detected source; as for the cumulative counts 
the horizontal error bars are not plotted as they are smaller than the size of the symbols. The Euclidean normalization is 
adopted so that for a static universe with a non-evolving population of objects and a constant luminosity function, the 
Euclidean normalized galaxy counts would follow a horizontal line. 

The differential number counts from the surveys at 15 /im (with fiux density converted into the 12 /im band, see 
previous Sect.) are also plotted in Pig. [3] (ELAIS-Sl, Pozzi et al. 2003; A2390, Altieri et al. 1999; HDP North, Aussel et 
al. 1999; HDP South, Marano PIRBACK Ultra Deep field, Marano Ultra Deep ROSAT field, Marano Deep field, Elbaz et 
al. 1999; Lockman Deep and Lockman Shallow field, Rodighiero et al. 2004). As for cumulative counts, the 15 /xm counts 
are corrected to 12 /im (see Sect. 3.1). Note that our faintest bin (at ~ 0.3 mJy) is very uncertain due to incompleteness, 
which may explain the downward shift for this last point. Por clarity we do not include the 12 /im galaxy counts from 
Clements et al. (1999) which show a large scatter, probably due to poor statistics (3-5 objects per bin). They may also 
suffer from contamination by a few stars at the bright end: the authors admit that there is at least one star in their 
galaxy counts (Clements et al. 2001). 

In Pig. [31 the ISO-ESS counts show a strong departure from Euclidean no-evolution models at faint fiuxes (< 1 mJy), 
with a very steep super-Euclidean slope. The same behavior is observed in all the other surveys plotted in the figure, 
showing agreement among the 15 /im surveys and with our 12 /xm survey. The low value of our last point at faint fluxes 
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Fig. 4. Example of the strong evolution in the SEDs predicted by PEGASE.3 for a star forming galaxy spiral Sc at 
various ages (increasing upwards), in the wavelength interval 1000 A to 1 cm. (see more details in Fioc et al. 2007) 

(~ 0.3 mjy) is due to the large uncertainty in the incompleteness correction. In the following, we address the issue of 
the origin of the excess counts by modeling the faint galaxy counts with the evolutionary code PEGASE. 

4. Modeling MIR galaxy counts per type with the code PEGASE.3 

In the following we adopt a flat Universe with the standard cosmological parameters: Hq = 72 km s^^ Mpc^^, Hm = 0.3, 
Ua = 0.7 (Spergel et al. 2003). 

4.1. The code PEGASE.3 

The new evolutionary code PEGASE.3 (Fioc et al. 2007) is an extension of the code PEGASE. 2 to the dust emission 
wavelength range (Fioc & Rocca-Volmerange 1997, 1999b, see also http://www2.iap.fr/pegase). The SEDs of reddened 
galaxies are consistently computed from the far-UV-optical-NIR to the FIR (far-infrared) domain. PEGASE.3 calculates, 
in a coherent manner, the stellar emission, extinction, metal-enrichment, dust mass and the emission of grains statistically 
heated by the radiation field. Two distinct dust media (interstellar medium and HII regions) are considered. As in 

PEGASE. 2, radiation transfer is computed in two geometries (slab and spheroid) appropriate for disk galaxies and 
ellipticals. Temperature fluctuations of the Polycyclic Aromatic Hydrocarbons (PAH) as well as graphite and silicate 
grain properties are derived with the method of Guhathakurta & Draine, (1989). For illustration. Fig. |3]shows the strong 
evolution of the spectral energy distribution (SED) of the spiral Sc template (specially in the MIR between 10 fim and 
20 fj,m) at various ages. 

4.2. Evolutionary scenarios of galaxies and relative number density fractions 
4.2.1. Model 1: "normal" evolved types 

In a first step, we adopt the same set of evolutionary scenarios of "normal" galaxies previously determined with the 
code PEGASE. 2 (Fioc & Rocca-Volmerange 1997) which fit the colors of nearby galaxies by types and the deepest 
multi- wavelength {Bj, U (and F30GW), / and K) faint galaxy counts in the UV-optical-NIR ranges (Fioc & Rocca- 
Volmerange 1999a). This set corresponds to the 8 following types: irregular magellanic Im; spirals Sd, Sc, Sbc, Sb, Sa; 
and ellipticals E2 and E). We use the same parameter set with the new version PEGASE.3 (Fig. [5]) and compute the 
continous SEDs of each type, extended to the mid- and far-IR, taking into account stellar and dust emissions as well as 
coherent absorption. The initial mass function is from Rana & Basu (1992) for each type. These templates are also able 
to predict photometric redshifts up to z—A with an accuracy az ~ 0.1 (Le Borgne & Rocca-Volmerange 2002). Therefore 
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the evolution scenarios are considered as robust. The main parameters (star formation law, initial mass function, galactic 
winds and astration rate) are listed for each type in Fioc & Rocca-Volmerange (1997) and Le Borgne & Rocca-Volmerange 
(2002). Star formation rates are proportional to the current gas mass density, a highly conservative assumption. The 
astration parameter varies with galaxy type. The current gas content -Mgas(0 is ruled by star formation, stellar 
ejecta, galactic winds and infall rates as described in Fioc & Rocca-Volmerange (1997); the adopted values by types are 
recalled in Table 1. 

In Table 1, we present the characteristic luminosities log(i*(12 /im)/^©) and log(L,(24/im)/LQ) adopted for the 
various galaxy types. We compute L^(12 iim.) and L*(24/xm) from the L^[Bj) values of the optical luminosity functions 
by types (Heyl et al. 1997), used to fit the faint optical counts (see Table 1 of Fioc & Rocca-Volmerange 1999a), and 
the colours Bj — 12 ^m and Bj — 24 /^m from PEGASE.3 at z = for each galaxy type (Fioc et al. 2007). The filter 
12 /im means /6'0/LWlO and 24 /im means Spitzer/MIPS 24 /xm; Bj is the blue Kodack Illa-J plus GG395 corrected filter 
(Couch & Newell 1980). The filter corrections from IRAS/12 ^m to ISO/LWIO and from IRAS/25 pm to Spitzer/2A fj.m 
are taken into account by the code. Among normal galaxies at z ~ 0, spirals Sbc are the brightest emitters at 12 /im 
and 24 /im, and also the most numerous (see Fig. [1}. We therefore assign to type Sbc the characteristic luminosities 
log(i,(12 /xm)/Lo) = 9.8 and log(i»(24/im)/io) — 9.6 derived from the observed z ~ luminosity functions measured 
by Rush et al. (1993) and by Shupe et al. (1998) respectively. We then scale the MIR luminosities of all the other types 
accordingly. 

Because the evolutionary scenario of ellipticals (see Fig. 3a of Rocca-Volmerange ct al. 2004) may play a specific role 
in the interpretation of observations of the ultra luminous infrared galaxies (ULIRGs), defined as galaxies with infra red 
luminosities >1O^^L0, it deserves more attention. The intense star formation rate (low p2 value) in the first Gyrs is fueled 
by a high infall rate from the gas reservoir. The activity is so intense that a huge dust mass is formed at early epochs 
from stellar ejecta, specifically from massive supernovae SNII. In normal elliptical galaxies, the star formation activity is 
supposed to be halted when strong galactic winds produced by the bulk of supernovae expel all the gas and dust content 
from galaxy. Most of stars and dust are already formed when the galaxy age is of about 1 Gyr, corresponding to z ~ 4 in 
the adopted cosmology. "Normal" ellipticals contribute very little thereafter to the infrared emission as they are largely 
dust-free from this age to the present time (0 < z <4). This scenario at z=0 matches the observation that the cold grain 
component (~ 50K) in elliptical galaxies has almost no contribution to the MIR flux (Xilouris et al. 2004). 

Column 7 of Table 1 gives the relative number fractions of galaxy types for model 1 , as derived from the UV-optical- 
nearlR. This model has no "ULIRG" component and is thus composed only of normal types (26.5% ellipticals (E -I- E2), 
23.7% Sa to Sbc spirals, 33.2% Sc, Sd spirals and 16.6% irregulars) which were found to fit the deepest UV-optical-nearlR 
faint galaxy counts (see Fig. 5 from Fioc & Rocca-Volmerange 1999a). 

4.2.2. Model 2: dusty ellipticals as ULIRGs 

Column 8 of Table 1 describes our model 2 with "ULIRG" which we use to adjust the MIR galaxy counts. In this model, 
1/3 of the ellipticals have over-luminous MIR luminosities given by log(L*(12 /zm)/^©) = 9.9 and log(L*(24 /zm)/^©) = 
9.7, which correspond to the observed characteristic L* of the MIR luminosity functions at these two wavelengths. At 
z = and in the MIR, they are nearly as bright as Sbc spirals, and respectively brighter by 2.5 mag at 12 pm and 5 
mag at 24 /im than normal ellipticals, whatever their type (E2, E). However, they better follow the evolution scenario 
of elliptical galaxies of type E, with the same astration rate P2^, infall and galactic winds at the same age. These 
overluminous ellipticals, forming at early epochs large masses of dust and stars, become much brigther at high z than 
spirals. This model does not need any additional starburst, as that seen in the typical case of M82 (Silva et al., 1998). 
Our model remains compatible with occasional starbursts of short time-duration (< 10® yrs), concerning a small mass 
fraction relative to the massive underlying elliptical galaxy. 

Column 8 of Table 1 lists the number fractions for model 2. The population of ULIRG/dusty massive ellipticals 
correspond to 9% of the total number of galaxies. The number of normal dust-poor ellipticals is then reduced to only 
17.5% (only 2/3 of the ellipticals observed in the visible). The rest of galaxies are normal; model 2 therefore respects the 
majority of fractions by type derived from the UV-optical-nearlR galaxy counts. 

4.2.3. The k{z) and e(z) corrections per type 

To calculate the apparent magnitudes at high z, the evolutionary e(z) and cosmological k(z) corrections are computed 
for each type, as in Rocca-Volmerange & Guiderdoni (1988), and are applied to the z = SEDs: 

m{z,t,) = Af(0,io) + Hz) + e{z) + 5 logio[i:iA(l + zf] + 25 (1) 

where Da is the angular diameter distance in Mpc, Da{1 + z)^ the luminosity distance, tz and to the cosmic times at z 
and z = 0; internal extinction is taken into account in our scenarios; no Galactic extinction term is added, because our 
deep survey is made in a region of weak galactic absorption. 

The same formation redshift of Zfor = 10 is arbitrarily adopted for all galaxy types, following the most distant galaxies 
discovered at z > 6 (Hu et al. 2002; Chary et al. 2005). 
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Type 


P2 


logio 


logio 


Colour 


Colour 


Model 1 


Model 2 




(Myr) 


[L*(12^im) 


[i.(24 /im) 


Bj - 12 fim 


Bj - 24 /im 


normal 


normal 






/Le] 


/Lq] 


AB 


AB 


only 


+ ULIRG 


ULIRG 


100 


9.9 


9.7 






0% 


9.0% 


E 


100 


8.9 


7.7 


-1.00 


-2.60 


9% 


0% 


E2 


300 


9.0 


7.8 


-1.20 


-2.82 


17.5% 


17.5% 


Sa 


1400 


9.0 


8.8 


-0.45 


+0.21 


7.9% 


7.9% 


Sb 


2500 


9.4 


9.2 


+0.37 


+ 1.28 


7.9% 


7.9% 


Sbc 


5714 


9.8 


9.6 


+1.38 


+2.30 


7.9% 


7.9% 


Sc 


10000 


9.6 


9.4 


+1.57 


+2.46 


16.6% 


16.6% 


Sd 


14286 


9.6 


9.4 


+1.61 


+2.48 


16.6% 


16.6% 


Im 


16000 


8.8 


8.5 


+1.57 


+2.42 


16.6% 


16.6% 



Table 1. Characteristic luminosity L■^, in units of L© at 12 /im (Column 3) and 24 /im (Column 4) of the luminosity 
functions used here, as a function of galaxy type (Column 1). Column 2 gives the value of the parameter p2, the rate 
defined by SFRi{t) — Mgas(t)/p2- Galactic winds occur at an age of 3 Gyr in ellipticals of type E and ULIRG, and at 1 
Gyr in ellipticals of type E2; there are no galactic winds in spirals. Infall time-scale of ULIRGs is 100 Myr as for ellipticals 
of type E or E2; it regularly increases for spirals from 2.8 Gyr (Sa) to 8.0 Gyr (Im) (see text for details and references). 
The L^,{^2 /im) for the normal types (other than ULIRG) are derived from the L^,{Bj) of the observed optical luminosity 
functions from Heyl et al. (1997) and the Bj — 12 /im and Bj — 24 /im colours, in the AB magnitude system, are computed 

at z ~ with PEGASE.3 (Columns 5 and 6). The value of L*(12/im) assigned to Sbc galaxies, the brightest emitters 
in this band, is that of the observed luminosity function measured with IRAS (Rush et al. 1993; Shupe et al. 1998); 
the same offset as for Sbc is then applied to all types. The last two columns show the adopted number density fraction 
by type: model 1 (normal galaxies only, column 7) uses the type distribution derived from the UV-optical-nearlR faint 
counts (Fioc et al. 1999a) while model 2 (normal galaxies + ULIRG, column 8) is built by replacing 1/3 of normal dust- 
poor ellipticals (9% of all galaxies) with ultra luminous elliptical galaxies (called ULIRGs). These last galaxies evolve as 
dusty ellipticals, they are ~ 2.5 to ~ 5 mag brighter than normal ellipticals (depending on wavelength) , and are thus as 
luminous as spirals Sbc at z = 0. The number densities for the other galaxy types are identical in both models. 



4.3. The 12fim luminosity function by type at z—0 

We take advantage of the quasi-similarity of the two filters IRAS /\2 fim and ISOCAM/JJWlO/12 /im (the flux differences 
are < 10%, see companion article) and use the local 12 /im luminosity function measured by Rush et al. (1993) from the 
IRAS catalogue (this measurement was later confirmed and extended to fainter fluxes by Fang et al. 1998), 

We adopt the values a — 1.7, (3 — 3.6 measured by Rush et al. (1993) for the non-Seyfert galaxies, which include all galaxy 
types detected at z = 0. As shown in Table 1, the i*(12 /im) adopted for each galaxy type are derived by combining the 
optical L^,{Bj) (Heyl et al. 1997) with the Bj-12 /im and Bj-2A /im colors predicted by the code PEGASE.3 (columns 5 
and 6). We assign to Sbc galaxies, the brightest emitters in these bands, the IRAS values of L*(12 /im) from the observed 
luminosity functions (Rush et al. 1993, Fang et al, 1998). The same offset as for Sbc is then applied to all types. 

5. Results of galaxy count modeling at 12 /xm 

5.1. Cumulative and difFerential counts from PEGASE.3 at 12/im 

Faint galaxy counts are computed following the already published formalism (see Sect. 2 in Guiderdoni & Rocca- 
Volmerange, 1990) which assumes that the number of galaxies of each type is conserved. Comparison of the ISO-ESS 
number counts with the predictions of PEGASE.3 at 12 /im is shown in Fig. [5] for cumulative counts, and in Fig. [6] for 
Euclidean-normalized differential counts (solid line in both graphs); the number counts from the 12 /im ISO-ESS, and 
from the published 15 /im surveys are also shown with the same line/symbol coding as in Figs. [5] and [3] Fig. [5]shows that 
the model with ULIRGs is in good agreement with the ISO-ESS cumulative counts. After colour correction, it also fits 
the deep ISO counts at 15 /im as well as the ultra-deep survey down to 0.05 mJy in the cluster-lens A2390 (Altieri et al. 
1999; Lemonon et al. 1998), implying that there is no significant number density variation between the field and clusters. 
The differential counts (Fig. [6]) , which are more constraining, remain in reasonably good agreement with the data, in 
particular at faint fluxes. The Euclidean-normalized differential counts predicted with PICGASE.3 do show the departure 
from the Euclidean cosmology (horizontal line) observed at 0.3 mJy in all the data samples. This bump is not due to the 
evolution of bright spirals, nor to normal early-type galaxies, but only to the evolution of the third of elliptical galaxies 
(9% of all galaxies) which are dusty ultra-bright ellipticals. Only from a few 0.1 to ~1 mJy, the slope of the Marano 
Deep Field is respected by models, in excess relative to other observations by a factor 2. The model prediction with only 
the universe expansion effect (obtained by applying only the k{z) corrections to the SEDs) is also plotted in Fig. [5] and 
Fig. [5] as a dashed line: it is noticeably insufhcient to reproduce the marked excess counts and the peak at 0.3 mJy. We 
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Fig. 5. The predictions of 12 /im cumulative galaxy counts (heavy line) calculated from model 1 of normal galaxies (left) 
and model 2 with 9% of ULIRGs (right): see Table 1, columns 7 and 8). They are compared to the observed 12 /xm 
cumulative counts from the ISO-ESS survey (black squares). The models with the cosmological k- correction only (no 
evolution correction: heavy dashed line) and the models with only the elemental comoving volume effect (dotted line) are 
also plotted for comparison. The faint counts from the 15 /im surveys are shown with the same line coding as in Fig. [21 




10-1 IQO IQl IQ2 10-1 10" 101 102 

12 fim fiux density {mJy) 12 /j.m flux density (mJy) 



Fig. 6. The predictions of Euclidean-normalized differential 12 ^m galaxy counts, calculated with model 1 (normal, left) 

and model 2 (with ULIRGs, right) using PEGASE.3 (solid hue), are compared to the observed differential counts at 12 
from the ISO-ESS survey (black squares). The cases without evolution (only cosmological k- correction) are shown with 
a dashed line, the case without fc + e corrections (only the comoving volume effect) with a dotted line. The symbol coding 
for the 15 /im surveys is the same as in Fig. [31 

also checked that the model without any correction (i.e. no k{z) + e{z) corrections applied to the SEDs), which only 
includes the evolution of the comoving elemental volume, yields decreasing differential counts which are incompatible 
with observations; this curve is shown in Figs. [HI and Fig. [HI as a dotted line. 

Note that the comparison of the observed ISO-ESS counts with the Euclidean case is more meaningful in the flux range 
where the number density is the highest. When galaxy numbers are statistically too small, error bars are large as shown 
at fluxes higher than 10 mJy and < 0.1 mJy. Finally we can ask whether the normal spirals which are ultra- luminous in 
the MIR (log(L* (12 /im)/L0) = 9.8 and log(i,(24/im)/L0) = 9.6 (see Table 1) can also reproduce the bump of MIR 
counts as well as models with ULIRGs. At 12 fim, the normal populations, dominated in the MIR by spirals Sbc are 
quite unable to reproduce the excess counts of the MIR surveys (Figs [SI and [51 left), even by adding 9% to the ^57% 
of normal spirals in place of gas-poor elliptical galaxies. 

6. Results of modeling Spitzer galaxy counts at 24 

6.1. The Spitzer surveys 

The deepest Spitzer/MIPS 24 /im surveys are in the area of the Chandra Deep Field South (Papovich et al. 2004). The 
corresponding galaxy counts (cumulative, shown in Fig. [71 and differential normalized to the Euclidean case, shown in 
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24 flux density (mJy) 24 /^m flux density (mJy) 

Fig. 7. Cumulative 24 /im faint galaxy counts (solid line) predicted by model 1 (normal, left) and model 2 (with ULIRGs, 

right) using PEGASE.3 (the same models as for the interpretation of the ISO 12 /im counts), are compared to the 
Spzfaer/MIPS/24 /im observations by Papovich et al. (2004). The luminosity function is from IRAS /2b fim (Shupe et 
al. 1998). As in previous figures, the dashed line is for counts taking into account only the k- correction and volume 
expansion effect, while the dotted line corresponds to the comoving volume correction only (no k + e correction). 




10-2 10-' 10" IQi 10-2 10-1 IQO IQl 

24 flux density (mJy) 24 /.im flux density (mJy) 



Fig. 8. The differential 24 nm faint galaxy counts computed with PEGASE.3 (solid line) for model 1 with normal galaxies 
only (left) and model 2 with ULIRGs (right), both compared to the SPITZER/MlPS/24: fim observations of Papovich et 
al. (2004). The other symbols are identical to the previous figure. For the 24 /im filter, the k-correction line (dashed- line, 
left) appears below the comoving volume line because the k-correction is negative, as shown by the SEDs displayed in 
Fig. 12 

Fig.[8|) are characterised by a typical bump of the galaxy density between 3 and 0.03 mJy, similar to that observed in the 
12 fim and 15 /im counts. This evolution signature of the 24 /im counts is confirmed by Marleau et al. (2004) and Chary 
et al. (2004). More recently, in the GOODS-ELAIS-Nl field, Rodighiero et al. (2006) have lowered the confusion limit by 
about 30-50% using a deblending technique, which leads to a decrease of the bright differential counts by a factor ~ 3, 
and an increase of the count slope at faint fluxes. The statistics of the 24 nm observations are poor for flux densities 
> 10 mJy and the samples suffer of incompleteness for fluxes < 80 /iJy. 

6.2. The local IRAS/25 iim luminosity function 

To model the Spitzer /2A fxra counts, we use the IRAS/25 //m luminosity function (corrected for 24 /j,m) as measured 
by Shupe et al. (1998), 

ML)=c(- + -^)x^-"{l + x)-'' , (3) 

\X 1 + X J 

where x = L/L^,, and the parameters are a = 0.437 and (3 — 1.749 for all galaxy types. As for 12 /im, we compute L*(24 
yum) by types from the Heyl et al. (1997) L^{Bj) and the Bj-2Aiim colour computed with PEGASE.3 (including filter 
correction from IRAS/25 /im to Spitzer /24: /im). We assign to spirals Sbc the value of L,(24 /im) derived from Shupe 
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et al. (1998). The same offset as for Sbc is applied to all types. The resulting values of the characteristic and the 
corresponding fractions per galaxy type are listed in Table 1. 

6.3. Cumulative and differential counts with PEGASE.3 at 24/im 

We model the faint galaxy counts through the Spitzer /24: /im filter with the code PEGASE.3 and the same evolving 
galaxy population (evolutionary scenarios, density fractions) as already used to predict the 12 /zm counts. Fig. [7] and 
Fig[5] present the comparison of the models 1 and 2 with the cumulative and differential number counts, respectively, 
obtained by Papovich et al. (2004). As for the 12 /xm counts, the 24 fj,m cumulative counts are well reproduced from the 
faintest flux up to a few mJy. Fig. [5] shows that the marked steepening of the differential counts normalized to Euclidean, 
and the subsequent decrease at faint fluxes are predicted by PEGASE.3, with a peak at ~ 0.3 mJy as observed. Note that 
the departure of the model from the observations at bright and faint fluxes is in good agreement with the data corrected 
for incompleteness and deblending (see Rodighiero et al. 2006). Once again, the fit is more meaningful in the flux range 
where the number density is the highest, as objects observed at high fluxes in the survey area are rare. Moreover, for the 
24 /xm filter, the fc-correction is negative at z < 2 as shown in Fig. |4]by the SED slope from 24 nm to ~ 8 /zm. In contrast 
with the study by Gruppioni et al. (2005), based on the flux density ratio S'24^m/'S'i5^m, our modeling does not require 
a population of additional starbursts, but rather the very strong evolution factor at high redshift of the star formation 
rate of elliptical galaxies as presented below. 

7. The cosmic star formation rate SFR(z) global per type 

The models of galaxy populations which simultaneously fit the 12 /im and 24 /im galaxy counts can be used to predict 
the cosmic star formation rate at high z. We then compute the global star formation rate SFR{z) as, 

SFR{z) = SFR,{z), (4) 

where SFRi(z) is the star formation rate SFRi{z) per galaxy type i. The total cosmic SFR{z) (heavy line) and SFR by 
groups of galaxy type SFRi(z) (thin lines) are presented in Fig. [5] in arbitray units (mass per year per volume unit). We 
distinguish the group 1. of E/SO galaxies which shows a striking SFR increase from z > 1 up to the redshift of formation 
(2^/or=10). The bulk of stars formed at early epochs (z > 4) while the SFR at z ~ is null because previous galactic 
winds have ejected all the gas content. Then the group 2. of early spirals (Sa-Sbc) which is dominant from z = 0.6 to 
z 4 with a predicted SFR which increases from z = to 1 by a factor of ~ 10, as found by Lilly et al. (1996), Madau 
et al. (1996), Connolly et al. (1997) (empty circles), mainly from the /-band selected CFRS observations. Finally the 
group 3. of late-type spirals and irregulars (Sc-Im) is dominant from z ~ to z 0.5 but has the faintest SFR at z > 2. 

Because ULIRGs follow the SFR scenario of ellipticals, both normal ellipticals and ULIRGs are represented by the 
same line (E/SO) on Fig. [9] This hypothesis, fully justified if the origin of the huge luminosity of ULIRGs is not stellar 
(i.e. if it is due to an active nucleus), emphasizes that the star formation history shown in Fig. [H] is a lower limit because 
it does not include the possible contribution from starbursts caused by interactions. Our modeling could accept some 
such occasional starbursts as long as they represent only a few percents of the total stellar mass and are short. 

The evolution of the SFR is different in the two groups Sa-Sbc and Sc-Im. Between z ~ 1.5 and z = 0, which 
corresponds to ages of ~ 6 to 13 Gyr, the star formation rate of a 10^^ M0 Sa spiral decreases from ~ 20 to 2 Mq yr~^ at 
z = 0. In the same redshift interval, the star formation rate of a less-evolved Sc spiral with the same mass increases from 
~ 3 to 4.5 M0 yr~^. This is consistent with the local observations: at z = 0, the current mean SFR of galaxies is ~ 5 
Mq yr~^ for late spirals of 10^^ Mq (and ~ 0.5 Mqut^^ for 10^" M0 irregulars); it is as low as 2 Mq yr^^ for early 10"'^^ 
Mq spirals (Kcnnicutt 1983). This SFR(t) variation is explained by the corresponding Mgas(t) variation. Because early 
spirals formed stars more efficiently in the past, their gas reservoir becomes depleted at z ~ 1.5 and the gas-dependent 
SFR is then rapidly decreasing by lack of fuel. This effect can explain the apparent "down-sizing" of galaxies (Panter et 
al. 2006). 

When summing over all galaxy types, the total cosmic SFR(z) (thick line) increases at high redshift: it evolves by 
a factor < 3 between z = and 1. This result is in agreement with the gradual decline of the UV luminosity from the 
deep multi- wavelength Hawaiian surveys (Cowie et al. 1999). This evolution is much smaller than the factor 10 decrease 
of the SFR derived from the less deep /-band selected CFRS sample (Lilly et al. 1996; Tresse et al. 2002). The CFRS is 
an /-band selected sample, thus biased towards early spirals, the major emitters in the / band. This is confirmed by the 
good match between the SFR{z) of the CFRS (empty circles) and that modelled for the Sa-Sbc galaxies in Fig. [9l We 
conclude that the /-band selection bias excludes late-type galaxies, too blue for detection in the /-band at the depth of 
CFRS. 

8. Discussion 

The faint galaxy counts at 12 fim in the ISO-ESS survey area show the same typical signature of evolution at 0.3 mJy as 
already found in the ISO/15 /im and Spitzer/2A /im galaxy counts. The careful flux calibration of the galaxy catalogue 
is based on the optical-MIR statistical properties of stars. Our results are consistent with the other MIR counts at 15 
/Ltm, thus demonstrating that the surveyed area of 800 arcminutes square is sufficient to average the inhomogeneities 
and properly analyze the populations of galaxies. 



Rocca-Volmerange, Lapparent, Seymour & Fioc: /50-ESS Counts at 12 



13 




Fig. 9. Histories of the cosmic star- formation rate density (global and by types) derived from models of galaxy populations 
which fit optical and IR faint galaxy counts. The total SFR(z) (heavy solid line), summed on all types is shown. The 
three SFRi{z) (thin lines) correspond to: 1) the major SFR at high reshifts {z >4) is from E/SO galaxies, due to the 
bulk of stars formed at early epochs, its contribution at z ~0 is null because previous galactic winds. 2) the intermediate 
SFR is from early spirals (Sa-Sbc), dominant from z = 0.6 to 4; it also fits the slope of the I-selected CFRS observations 
(empty circles), up to z=l before incompletness. 3) the weakest SFR are from Sc-Im galaxy types (late- type spirals and 
irregulars) except at z ~0 where they are numerous and actively forming stars. Models of ULIRGs follow the scenario of 
ellipticals, strongly evolving from z >2 to Zfor (here arbitrarily adopted at 10). This is the most conservative hypothesis 
for ULIRGs: ultra-brightness is likely attributed to active nucleus while adding occasional starbursts would increase the 
global SFR{z). 



The most important result of our analysis is that a minor (< 10%) population of dusty ultra-bright elliptical galaxies 
can explain the excess of the mid-IR emission observed in the 12 /xm and 24 /im faint galaxy counts at ^ 0.3 mJy. Here, 
due to its high IR brightness, this population is associated to ULIRGs, while the other populations, seen in the UV- 
optical-nearlR as well as in MIR counts, are called normal galaxies. Because the evolutionary code PEGASE.3 predicts 
multi- wavelength SEDs by simultaneously following the evolution of stars, gas, dust and metal-enrichment, our analysis 
results find a natural explanation in the basic scenarios of galaxy evolution. 

The strong advantage of our analysis is the multi-wavelength approach: the same SFR scenarios already found to fit 
the UV-optical-nearlR galaxy counts (Fioc & Rocca-Volmerange, 1999a) are also applied to the MIR (12 /im, 15 /im 
and 24 /im) using ISO (see companion article and Elbaz et al. 1999) and Spitzer satellites (Papovich et al. 2005; Le 
Floc'h et al. 2005). These evolutionary scenarios are robust because the evolution time-scales of the dominant emitters 
at various wavelengths (massive stars, evolved stars, dust grains from the interstellar medium and HII regions) go from 
a few million years to ~ 13 Gyr. In the mid-infrared, the new model with ULIRGs, proposed to fit ISO/12 /im galaxy 
counts does not likely change the UV-optical-nearlR predictions; moreover it is confirmed by the ISO/lbfiia and the 
more recent Spitzer /2A pm galaxy counts. 

One difficulty of the interpretation is that, at z ~ 0, the brightest Sbc spirals appear as bright as dusty ellipticals in 
the MIR. Fig. [5] to Fig. [5] show that the populations of "normal" ellipticals, spirals and irregulars (model 1) seen in the 
optical are largely insufficient to reproduce the 12 /im and 24 /im differential and cumulative counts. Even the brightest 
spirals (with log(L*(12 yum)/L0) = 9.6), which have IR luminosities comparable to ULIRGs at z = 0, decline too rapidly 
at increasing z to fit the 12 fun and 24 /im number counts. 
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The surprising result is that we succeed in reproducing the typical excess of MIR counts observed at c± 0.3 mJy by 
only replacing 9% of "normal" galaxies in the optical (and 1/3 of the ellipticals) with ultra- bright galaxies in the IR. In 
fact at high redshift, the evolution correction takes over as the major parameter: it is noticeably insufficient for spirals; 
only ellipticals have sufficient star formation rates to reproduce the stellar and dust emission at high redshifts. The 
fraction of these objects is small and the enormous luminosities required can only be reached if these ULIRGs contain 
huge dust masses heated by large numbers of energetic photons. 

The star formation history of ULIRGs follows that of elliptical galaxies in Fig. [9l large masses of stars and dust are 
formed at high z. ULIRGs would appear as "normal" ellipticals in the optical with masses of ~ 10^^ M©; they could 
even be more massive if the excess of stellar luminosity is hidden by an excess of dust. 

This dusty massive population, revealed at high z in the mid-infrared, evokes the population of high-z (> 4) massive 
ellipticals found in the K-z Hubble diagram (Rocca-Volmerange et al. 2004), and confirmed in the rest-frame H-z diagram 
with Spitzer (Seymour et al., submitted): these distant galaxies are forming also high masses of stars and dust. Due to 
the short time-scales required to build such objects at z > 4 (< IGyr), both populations likely formed at early epochs 
by rapid dissipative collapse or major merging, rather than by slower hierarchical merging which would take ~ 10 Gyr. 

This interpretation of ULIRGs as dusty massive ellipticals is in agreement with the drastic evolution of the infrared 
luminosity function when compared to the UV luminosity function (Takeuchi et al. 2006). It is also in agreement with 
the result that the MIR-selected galaxies contribute to more than 70% of the Cosmic Extragalactic Background in the 
MIR (Dole et al. 2006): our modelling confirms that other galaxy types are too faint at both 12 /xm and 24 /xm. However, 
the conclusion of the last authors that galaxies contributing the most to the total cosmic infrared background have 
intermediate stellar masses is not confirmed by our results. 

One puzzling issue is to know how the dust mass could be maintained inside the galaxy host, hence no released by 
galactic winds according in normal ellipticals. To maintain the large dust mass of ULIRGs/dusty ellipticals within galaxy, 
dust must not be mixed with gas and stars. It has to be located in preferential zones such as the galaxy core. In the core, 
where the potential well becomes intense in case of an embedded black hole, dust could be preferentially retained. In the 
AGN environment, the deep potential well would drag dust more efficiently, as it is more massive than gas; dust then 
would fall more rapidly down in the inner core. Moreover if a Compton thick AGN is embedded within the proposed 
ULIRG/dusty ellipticals, the large variability from one ULIRG to another (Armus et al. 2006) is explained by orientation 
effects. Only observations at high spatial resolution will allow the dust geometry to be determined. 

The other issue is the presence of a large number of energetic photons heating dust grains at all ages. To produce 
them, massive stars and/or the presence of an AGN can be evoked. Our results do not exclude that, in the case of galaxy 
interactions, an exceptionally extincted starburst, undetected in the optical, could be ultra-bright in the IR. But such 
an event is rare and is not representative of a galaxy population on a long time-scale. Note that despite the high value 
of their MIR luminosity (log(L*(12 /im)/LQ) = 9.9) at z = 0, the proposed ULIRGs faintly contribute to the faint UV- 
optical-NIR galaxy counts by their number. They also are suffering an exceptional extinction due to their dust amount. 
Further spectral syntheses and high spatial resolution are clearly needed. The recent analysis by Takeuchi et al. (2006), 
based on the combination of data from the UV satellite GALEX and from IRAS shows that the luminosity function 
evolves more strongly in the far-infrared than in the far-UV. This is compatible with our dusty elliptical population. At 
last, by analyzing the far-UV galaxy counts with FOCA at 2000A, Fioc & Rocca-Volmerange (1999a) suggested that a 
fraction of episodic starbursts could be required to interpret the UV excess of galaxy counts, in addition to the normal 
populations of galaxies. However, the number density, weak star formation rate, and low metallicity of these populations 
are not sufficient to explain the excess of MIR luminosity at high redshifts. 

Finally, these scenarios of ultra-luminous galaxies at high redshift, imply a very rapid phase of mass accumulation. This 
is also supported by the fact that ULIRGs, evolving as ellipticals and hosting a hidden AGN, look like the population 
of radio-galaxies discovered from the K-z diagram at high redshifts which also show strong and hot dust signatures 
(Rocca-Volmerange & Remazeilles 2005). However, the proposed population of ULIRGs derived from the infrared is 
more numerous (9%) than the radio-galaxy hosts detected in the optical (< 4%). This indicates that half of the ULIRGs 
could be so obscured in the optical that they would be invisible. They may however be revealed at other wavelengths. 
Several surveys have discovered populations of AGN which appear brighter and denser than the classical populations 
identified in the optical (see for example Martini et al. 2006). The population of hyper-LIRGs (L//j > lO^^L©), sometimes 
associated to Ly^ blobs, has a low LLy^/iboi efficiency (0.05-0.2%) according to Colbert et al (2005). The 12 /^m and 
24 /im galaxy counts analysed here may correspond to the best wavelength domain where this population of possibly 
embedded AGNs could be detected. 



9. Conclusion 

We present the faint galaxy counts derived at 12 /im from the observation of a large and deep mid-infrared (MIR) survey 
in the field of the optical ESO-Sculptor Survey (ESS), through the large LWIO filter with the ISOCAM instrument on 
board the ISO satellite. The infrared observations cover an area of ~ 75% of the ESS spectroscopic survey, where galactic 
cirrus is sparse, and were performed in continuous raster mode. The fiux calibration has been adjusted by using the 
optical-infrared IRAS colours {B — 12 /im, V — 12 /tm) of standard stars. Because of its large area of ~ 680 arcmin^, 
the /S'O-ESS survey provides complete 12 /tm galaxy counts down to 0.24 mJy, after incompleteness corrections (using 
two independent methods). The full data analysis and the resulting catalogue of 142 detected sources is published in the 
companion paper, Seymour et al. (2007). 
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The galaxy counts are presented using two different binnings: cumulative counts iV(> S), to avoid the fluctuations in 
the number density per bin, and Euclidean normalized diS'erential counts -S'^-^dN/dS. When corrected for incompleteness, 
both the cumulative and differential ISO-ESS 12 /xm counts averaged over the ~ 680 arcmin^ area show good agreement 
with the existing measurements in the closc-by ISOC'AM filter at 15 /im, after correction for the different wavelengths. 
In particular, the Euclidean-normalized differential counts of the ISO-ESS survey display the same excess as the other 
existing MIR surveys at flux densities of 0.3 mJy. This excess is also observed in the Spitzer galaxy counts at 24 /um. 

We propose an interpretation of the cumulative and differential counts with the help of the new evolutionary code 
PEGASE.3 (Fioc ct al. 2007). For each galaxy type, PEGASE.3 predicts the spectral energy distributions from the 
optical to the far-IR; the emission of stars and dust, the extinction, star formation history, metal enrichment and dust 
mass are computed consistently. 

With these evolutionary standard scenarios we have successfully modelled the multi-wavelength faint galaxy counts 
in the far-UV, optical and near-infrared (Fioc & Rocca-Volmcrange 1999a). In the present article, we are able to simul- 
taneaously flt the ISO-ESS 12 /im, ISO 15 /im and the Spitzer2A /im faint counts, by increasing by ~ 2.5 mag (at 12 pm) 
to 5 mag (at 24 /im) the luminosity of a small fraction of galaxies (9%; all of elliptical type), while the rest of galaxies 
(17.5% normal ellipticals, 57% spirals and 16.5% irregulars) are identical to the galaxy populations already known from 
the UV-optical-NIR surveys. The ultra-bright galaxies display all the characteristics of ULIRGs and appear as "normal" 
ellipticals in the optical. Because these results cover a very large wavelength domain, from the UV-optical-NIR to the 
MIR (12 /xm, 15 pm and 24 /im), they confirm the robustness of our scenarios. 

The other important point is that no additional starbursts are required to fit the MIR excess. Highly luminous 
starbursts with a short e- folding time (~ 10'' years), which may explain some nearby ULIRGs such as M82, are not 
incompatible with our results if they remain exceptional objects. Another point is that the normal galaxy populations, 
including the bright IR spirals, are insufficient at high redshifts to fit either the cumulative or the differential ISO-ESS 
12 fim and SPITZER 24 /xm counts. 

The star formation history of the proposed ULIRGs and normal ellipticals, in respective proportions 1/3 and 2/3, 
can fully explain the excess in the cumulative and differential galaxy counts at 12 /xm and 24 /xm, thanks to a large dust 
mass produced by early-formed stars. We suggest that in these ULIRGs, dust is not mixed with stars, as otherwise it 
would be expelled by the galactic winds. The possibility of an embedded Compton thick AGN, which would explain the 
ultra-brigthness by keeping dust within galaxy, will require confirmation with other observations. Similar to the distant 
radio galaxies found in the K — z diagram (Rocca-Volmerange et al. 2004), the most distant ULIRGs would appear as 
the most massive ellipticals. As concluded in the mentioned article, this new result in the MIR favors the hypothesis 
of a galaxy evolution process based on dissipative collapse or a rapid merging on a short time-scale (< 1 Gyr) at high 
redshift (z > 5). 

We emphasize that these results are robust since they are derived from the simultaneous adjustment of the MIR 
12 /tm and 24 /tm counts respecting the majority (> 90 %) of galaxies observed in the optical faint number counts (down 
to B > 29 from the HDF-N, William et al. 1996). Moreover, our analysis by type allows us to identify the various factors 
explaining the steep increase of the faint galaxy counts. The comoving volume element effect, as well as the expansion 
effect (fc-correction), arc insufficient to explain the faint count excess. Higher spectral and spatial resolution observations 
associated to deeper counts at longer wavelengths with the Spitzer satellite and the future Herschel satellite will hopefully 
allow the detection of the dust emission emitted by early elliptical galaxies in their primitive epochs, validating our present 
results and allowing to hopefully observe primeval galaxy populations down to the deepest extragalactic backgrounds. 
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